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Abstract. The Heun equation can be rewritten as an eigenvalue equation for an ordinary 
differential operator of the form —d'^/dx^ + V{g] x), where the potential is an elliptic func- 
tion depending on a coupling vector g G K^. Alternatively, this operator arises from the 
BCi specialization of the BCn elliptic nonrelativistic Calogero-Moser system (a.k.a. the 
Inozemtsev system). Under suitable restrictions on the elliptic periods and on g, we asso- 
ciate to this operator a self-adjoint operator H{g) on the Hilbert space Ti = i^([0, wi], dx), 
where 2a;i is the real period oiV{g\x). For this association and a further analysis of H{g), 
a certain Hilbert-Schmidt operator I{g) on Ti, plays a critical role. In particular, using the 
intimate relation of H{g) and X(g), we obtain a remarkable spectral invariance: In terms of 
a coupling vector c £ R** that depends linearly on 9, the spectrum of H(g{c)) is invariant 
under arbitrary permutations a{c), a E S4. 
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1 Introduction 

In the first two papers of this series ([HE], henceforth referred to as I and II), we initiated 
a detailed study of spectral properties of operators arising in the context of elliptic Calogero- 
Moser type systems. The operators at issue are not only the Hamiltonians defining these systems, 
but also families of Hilbert-Schmidt (HS) operators that seem unrelated at first sight. The 
relation consists in the kernels of the HS operators being eigenfunctions of differences of the 
Hamiltonians. As we have first detailed in our lecture notes [3] and also sketched in I, this 
relation can be exploited to derive novel information about the Hamiltonians, viewed as Hilbert 
space operators (as opposed to differential or difference operators). 
The present paper is concerned with the Heun Hamiltonian 



H{g;x) 



+ y{9;x), g 



,91,92,93) G K^, 



V{9;x) = "^gtigt - l)p{x + uJt;'K/2r,ia/2), r, a > 0, 



(1.1) 
(1.2) 



t=o 



*This paper is a contribution to the Proceedings of the Workshop "Elliptic Integrable Systems, Isomonodromy 
Problems, and Hypergeometric Functions" (July 21-25, 2008, MPIM, Bonn, Germany). The full collection is 
available at |http: / /www. emis.de/journals / SIGMA /EUiptic-Integrable-Systems.html , 
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where 

loq = 0, cJi = 7r/2r, 1^2 = ia/2, LJ3 = ~loi — 102- 

(For a lucid account of elliptic and allied functions we recommend [1], whose conventions and 
notation we mostly follow.) It can also be viewed as the BCi 'nonrelativistic' elliptic Calogero- 
Moser Hamiltonian, the BCn case being the Inozemtsev system [5]. There exists an exten- 
sive lore about the eigenfunctions of H{g) viewed as a differential operator, cf. e.g. the re- 
views [UlTllH]. In recent years, this lore has been significantly enlarged by Takemura, whose 
work can be traced from his recent paper [S]. 

Here, however, we are concerned with spectral properties of H(g), reinterpreted as a self- 
adjoint operator on the Hilbert space 

n = L^{[0,TT/2r],dx). (1.3) 

A few of the properties we derive (such as the occurrence of solely discrete spectrum for suitably 
restricted g) follow from well-known arguments in the theory of Schrodinger operators [10] (or, 
alternatively, from Sturm-Liouville theory and the theory of eigenfunction expansions). Our 
main tool for reobtaining known features and obtaining novel ones is the integral operator 
whose kernel is given by 

■f{g; X, y) = w{g; xf'^S{g- x, y)w{g' ■ yf'^, x,y e (0, 7r/2r). (1.4) 

This kernel and the definition of its ingredients can be found in I, cf. I Subsection 3.2. For 
ease of reference, however, we repeat the relevant definitions: 



g' = JnQ, Jn = t, 



(I 1-1 l\ 
111-1 
-1111 

V 1 -1 1 1 / 



(1.5) 



Sig;x,y) = exp{-Sg\n[R{x + y)R{x - y)]), (1.6) 



^9 

^3 ,3 



t=o t=o 

w{g; x) = l/c{g; x)c{g; -x), (1.8) 

c{g; x) = R{x + ia/2)-3''R{x + ia/2 - -k /2r)-3^ R{x)-3^ R{x - ■K/2r)-3\ 
00 

R{z) = YI{1 - exp[-(2/ + l)ra + 2irz]){l - exp[-(2/ + l)ra - 2irz\), z G C. 
1=0 

At this point we would like to mention that a close relative of the kernel S{g\ x, y) has appeared 
in the literature before, a fact we only recently became aware of. Indeed, it can be tied in with 
the kernel that can be found on p. 124 of a monograph by Slavyanov and Lay [llj . (It seems 
this kernel was first introduced in a paper by Kazakov and Slavyanov [12].) More precisely, 
when a p-function substitution is made in the latter kernel, it can be factorized in terms of the 
Weierstrass cj-function; the nontrivial part is then an analytic continuation of S{g;x,y). Both 
in [11] and in Takemura's paper [13], this kernel is used to relate Heun functions defined on 
different intervals. 

The entire function R{z) plays a crucial role in this paper. It satisfies 

R{z) = ew(-f:^^^^^^], \^z\<a/2, (1.9) 
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R{2z)= Yl R{z + i5a/A)R{z + Ti/2r + i5a/A), (1.10) 
&=+- 

oo 

R{z + ia/2)=ipe-''''s{z), p = 2r\\{l - e-'^^'^'^f . (1.11) 

k=l 

Here, s{z) is related to the (T-function by 

s{z) = exp (— T^z^r/vr) (7(z; 7r/2r, za/2), (1-12) 
so that 

f{z + ia) = — exp(ra — 2irz)f{z), f = s, R. 

(For more details, see |14j . where the functions s and R were introduced.) 
From (II. lip it follows that the weight function (II. 8p can also be written 

w{g;x) =p'^3'^+^3is^xf90s^^/2r-xf9iR(xf9^R(-jT/2r-xf33_ (1.13) 

It is positive on (0, 7r/2r) for all g € M^, but integrable on [0,7r/2r] only for go,gi > —1/2. 
It now readily follows that for ^{g;x,y) to belong to the kernel Hilbert space 

nK = LmO,7r/2rf,dxdy), (1.14) 
the necessary and sufficient condition is that g belong to the parameter set 

Hee {^EM^ I go,gi,g'o,g[ > -1/2}. (1.15) 
Hence it is clear that the integral operator 

(J(g)/)(x) ^ r^^' dy^ig; x, y)f{y), f G Ti, g G U, (1.16) 
Jo 

is well defined and HS. 

It is not obvious, but true that its kernel is a zero-eigenvalue eigenfunction of a difference of 
Heun Hamiltonians: 

{H{g;x)-H{g';y))^{g;x,y)=0, x, y G (0, 7r/2r). (1.17) 

(Indeed, this is a corollary of Proposition 3.4 in I, cf. 1(3.115).) At face value, this eigenfunction 
identity seems to have no bearing on the Hilbert space eigenfunctions of the two (generically) 
distinct Hamiltonians. A principal result of this paper is, however, that for g in the restricted 
parameter set 

Ur = {geU\sg> 0}, (1.18) 

the Hilbert space eigenvector ONBs (orthonormal bases) of H{g) and H{g') are given by func- 
tions en{g; x) and en{g'; x), resp., which yield the singular value decomposition of I^g) (cf. [15]): 

00 00 
I{g) = ^M9)en{9)^en{g'), z^o > 1^1 > • • • > 0, ^z/2<oo, (1.19) 

n=0 n=0 
{em,en) = Smn, m, 71 G N. 

(All singular values are positive, since X(g() has trivial kernel and dense range for g £ H,., cf. ()2.4p . 
To appreciate the restriction on Sg, note that I{g) has rank 1 for Sg = 0, cf. (jl.6p .) 
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This result entails in particular that H{g) and H{g') have solely discrete spectrum, the two 
spectra being equal. To be sure, in the case at hand it follows from well-known Schrodinger 
operator theory that H{g) can be reinterpreted as a self-adjoint operator that has an ONB of 
eigenvectors, provided g is suitably restricted. The main novelty is that these eigenvectors are 
just those occurring for T{g), and that this fixes the self-adjoint extensions of the Schrodinger 
operator H{g) in case it is not essentially self-adjoint on Co°((0, 7r/2r)) (namely, for g^ and/or gi 
in (—1/2, 3/2)). As a bonus, this yields a hidden spectral invariance under permutations of four 
new parameters, which implies spectral consequences for special cases that are quite startling. 
Last but not least, the state of affairs for the present case gives rise to a simple paradigm for 
the 'relativistic' (analytic difference operator) Heun case, where the usual self-adjoint extension 
theory is of no help and no previous spectral results are known, cf. the fourth paper in this 
series [16]. 

We proceed to sketch the organization and results of this paper in more detail. We begin 
Section [2] by deriving features of functions in the range of T{g) , assuming (7 G 11^ . Then we 
show that I{g),g G 11^, has dense range and trivial kernel. These crucial features readily follow 
from Lemma l2.ll which asserts that the integral operator with kernel S{g;x,y) has kernel {0} 
on L^{[0,TT/2r]) whenever Sg is positive. The characteristics of the range of Tig) are now used 
as a guide for a choice of domain for H{g), namely the dense subspace D{gQ,gi) given by ()2.7p . 
Lemma |2. 21 then asserts that H{g) is symmetric on D{gQ,gi) for all g in the parameter space 

U = {gGR^\go,gi>-l/2}. (1.20) 

The proofs of these lemmas involve technicalities that are not enlightening, so we have relegated 
them to Appendix lAl 

Next, using the key identity (jl.lTp . we show in Lemma 12.31 that for all g G 11 the subspace 
D{go,gi) is a core (domain of essential self-adjointness [10]). This prepares us for studying the 
self-adjoint closure (again denoted H{g)) by using ^^{g). Only a few simple steps are then needed 
to reach the conclusion that the above ONB {6^(5)} may be chosen so that it is an eigenvector 
ONB for H{g). Moreover, since the eigenvectors belong to D{gQ,gi), it easily follows from the 
Frobenius method that the spectrum of H{g) is simple, cf. Theorem 12.41 

Whenever the potential V{g; x) is bounded below, it is obvious that H{g) is also bounded 
below. For gQ and/or gi in (0, 1), however, V{g;x) is not bounded below, so an extra argument 
is needed to show that for the pertinent self-adjoint extensions of (II. ip semi-boundedness is 
preserved. We prove this by a comparison argument involving the trigonometric BCi Calogero- 
Moser Hamiltonian 

(f 

Ht{go,gi;x) = -— ^ + Vt{go,gi;x), 
ax'^ 

\n ^ _ ^^go(go - 1) , r^g\{9\ - 1) 

yt\9Q,9\\x) = \ 2 • (1-21) 

sm rx cos^ rx 

The point is that this operator is well defined and essentially self-adjoint on D{go,gi), and 
that it has an ONB of eigenvectors in D{gQ,gi) involving the well-known Jacobi polynomials, 
cf. Theorem 12. 5[ 

We would like to add that the trigonometric Hamiltonian p.2ip plays a key role in Takemura's 
paper [T7], where he studies the Heun Hamiltonian H{g;x) (jl.ip with gQ,gi > 1 from a Hilbert 
space point of view. His strategy is based on a perturbation expansion in the parameter p = 
exp(— or). (In this connection it should be pointed out that H{g;x) with 52 = 53 = reduces 
to Ht{go,gi;x) for a 00.) 

The main result of Theorem 12.61 is that for g £ Ilr the spectra of H{g) and H{g') coincide. 
Unfortunately, we have not been able to show that the singular values in (|1.19p are distinct for 
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all g G Ilr, a feature that seems plausible to us. If this is indeed the case, then it would follow 
that en{g) is the i/(5)-eigenvector with eigenvalue En{g) in the obvious eigenvalue list 



In Section [3] we fix g G 11^. and address the question for which ^ G 11 the spectrum of H{g) 
coincides with that of H{g). From Theorem 12.61 we know this is true for g = g' , but there are 
a few more involutory (^-transformations for which this is valid. We determine the group G 
generated by these involutions in Theorem 13.11 It leaves the parameter set 



invariant, and is isomorphic to the permutation group 54. (For ease of reference, we collect the 
four parameter sets ([IM]), P^TS]) . pTT5]) and (fOO]) in a chain 

Hg c n,. c n c n, 

and note that each inclusion is proper.) Indeed, defining new couplings 

Co = 50 + 53 - 1, ci= gi + g2-l, £2 = 51-52, C3 = go-g3, 
the spectrum of 



is invariant under arbitrary permutations of cq, ci, C2, C3. 

In SectionlHwe consider various special ^-choices. For the sixteen (7- values yielding V{g; x)=Q, 
all of the relevant quantities can be determined in great detail. A striking consequence of our 
spectral invariance results is that there exist quite a few (7-choices for which V{g; x) 7^ 0, yet the 
spectrum of H{g) equals that for one of the cases where V{g; x) = 0. We also consider the case 
Sg = 0, for which T{g) has rank 1. In particular, using results from I and [18], we verify that I{g) 
connects the ground state of H[g') with the one of H[g). Finally, we add some observations on 
the case where the 7r/r-periodic potential V{g]x) has no poles for real x, a well-known setting 
that is for instance studied in considerable detail in |19j . 

The paper is concluded with Appendix [XJ in which we present the proofs of Lemmas 12.11 
and [121 

2 The eigenvector ONBs of H{g) and H{g') 

Unless explicitly mentioned otherwise, we assume from now on that the coupling vector g be- 
longs to the restricted parameter set 11^ (jl.lSp . This entails in particular that the integral 
operator T{g) (jl.l6p is a well-defined HS operator on the Hilbert space TL (jl.Sp . We begin this 
section by obtaining features of functions in the range of T{g) . 
Consider first for f £ TC and x G [0, 7r/2r] the function 



Eo{g) < Ei{g) < E2{g) < ■ ■ ■ 



5 e n^- 



Hg = {5 G I 50 + 51 + 2 > 52 + 53} 



(1.22) 



H{c) ^ H{g{c)) 




^w{g';yf'^f{y) 



(2.1) 




""^(5'; 0^^^ £ deduce from the Schwarz inequality 
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It readily follows that h{x) extends to a function h{z), 5Rz = x, that is analytic in the strip 
\Qz\ < a/2. Also, h{z) is vr/r-periodic and even. Now R{z) is an entire, vr/r-periodic and even 
function without zeros for \Qz\ < a/2, so the function 

s(x)~^^«s(7r/2r-j;)~fi(j(g)j)(2;) =i?(a;)f2i?(^/2r-j;)53/i(x), xG (0,7r/2r), (2.2) 

extends to a function that is analytic for {'^zl < a/2, vr/r-periodic and even. 
Next, we point out the relations 

geUr ^ g' e Ur, 

I{gr=I{g'), (2.3) 

which are easily verified. As a consequence, the orthocomplement of the range of I{g) equals 
the kernel oiI{g'). From (j2.ip it is immediate that X{g),g G 11, has rank 1 when Sg vanishes. 
By contrast, for > the integral operator T{g) has dense range and trivial kernel: 

Ker(T((7)) = {0}, Ran(J(5))- = H, Vr? G H,. (2.4) 

This is clear from the above and the following lemma, which is proved in Appendix [Al 

Lemma 2.1. Letting Sg > 0, suppose (j){y) £ L^([0, 7r/2r], dy) satisfies 

fn/2r 

dyS{g;x,y)(j){y) = 0. 



Then (j) = 0. 

To proceed, we note that the kernel function ^{g;x,y) is real- valued. Prom (j2.4p . (j2.3p and 
the HS property we then deduce that for a fixed g G 11^, there exist two ONBs {en{g')}n&i and 
{en{g)}n&i of real- valued functions such that 

oo 

'^{9',x,y) =^Vn{g)en{g;x)en{g']y), i^o > z^i > • • • > 0, (2.5) 

n=0 

with {fn} G ^^(N). (The choice of the first ONE uniquely determines the second one; later 
on, we will return to the specific choice we make.) Moreover, since en{g) belongs to the range 
of T{g) , the above analysis implies that it is of the form 

en (5; x) = w{g; xf^'^dn{g; x), (2.6) 

where dn{g;x) extends to a function that is analytic in \Qz\ < a/2, even and 7r/r-periodic. 

We are now prepared to deal with the Heun Hamiltonian H{g). We choose as initial domain 
for H{g) the space 

D{go,gi) = s{xyos{Tr/2r-x)9'V, go,gi>-l/2, (2.7) 

where 

V^{feC\[0, 7r/2r]) | /'(O) = /'(vr/2r) = O} . (2.8) 

This domain choice ensures that the ONB functions en{g', x) belong to the domain. Indeed, from 
the paragraph containing (|2.2p we see that the function on the r.h.s. of (j2.2p belongs to T>, so 
that we have more generally 

RaniI{g))cD{go,gi). (2.9) 
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From en{g) G D{gQ,gi) it already follows that D{go,gi) is dense. But this feature can also be 
derived quite easily without invoking 1(g)- Indeed, T> contains the functions cos(2nrx), n G N, 
which are well known to yield an orthogonal base of 7i. Therefore, the assumption / _L D{gQ,gi) 
entails 

7r/2r 

dxs{xy°s{Tr/2r - x)^^ cos(2nrx)/(x) =0, Vn G N. 

From the orthogonal base property we now deduce f{x) = 0, so that D{go,gi) is dense. 
We are now prepared for our second lemma, whose proof is relegated to Appendix |X1 

Lemma 2.2. Assume g belongs to H ()1.20p . Then H[g) (jl.ip is well defined and symmetric on 
the domain D{gQ,gi) given by ()2.7p . ()2.8p . 

The proof of this lemma does not involve the HS operator 2{g)- In the next lemma we do 
use I{g) to obtain a short proof of essential self-adjointness on D{gQ,gi). 

Lemma 2.3. Assuming 5 G 11, the domain D{gQ,gi) is a core for H{g). 

Proof. We begin by noting that the potentials p{x + ia/2) and p(x — 7r/2r — ia/2) yield 
bounded self-adjoint operators, so that they do not influence domain issues. Therefore, we need 
only prove the core property for one value of 32 and 5(3. We choose 52 = 53 = 1 (say), the point 
being that this entails g G 11^. (This is easily checked, cf. ()1.5p . (|1.15p and (jl.lSp .) Hence we 
can now use T{g). 

Assume (p € TC satisfies 

{cl),{H{g)+i)i;)=0, VVgD(5o,5i)- 
By virtue of ()2.9p . this implies in particular 

{(t>,{H{g)+i)l{g)f) = 0, yfen. 
Thanks to (I1.17|) and (12. 3p . this inner product can be rewritten as 

{{H{g')-i)i{g')<pj) = o, yfen. 

This implies that the vector I(g')(j), which belongs to D{gQ,g^) by (12. 9p . is an eigenvector with 
eigenvalue i for the symmetric operator H{g'). Thus we must have I{g')(j) = 0. By (j2.4p . this 
entails (p = 0. Therefore Ran {H{g) + i) is dense in 7i. Likewise, Ran {H{g) — i) is dense in 7i. 
Hence the core property results. ■ 

We would like to point out that the dense subspace 

C = C7o~((0,a;i)) cZ)(<7o,9i), (2.10) 

is already a core for H[g) whenever go^gi > 3/2. This is basically known, but for completeness 
we sketch a proof. Choosing at first g G M^, consider the differential operator H{g) restricted 
to C. This is obviously a well-defined symmetric operator on H, and it is easy to verify that the 
domain of its closure contains the space 

C' = {i;e CH[0,iOi]) I ^(0) = iPiiOi) = V'(0) = V'(a;i) = 0, 

i;'{x) G C\{0,LJi)),i;"{x),x-^^{x),{uJi-x)-^ij{x) G Ti}. 



Now assume (p is orthogonal to [H{g)+i)C. Then (/> is a weak solution to the ODE [H{g)—i)(j) = 
on (0,a;i), and by hypo-ellipticity it is a classical solution. 
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Choosing next go,gi > 3/2, it follows from the Frobenius theory that near we have 
(pix) = x^«/io(x), 
with /lo(a^) analytic at x = 0. (Indeed, 

2,1-90 ig 

not in L near 0.) Likewise, near x = cJi, 

<Pix) = (a;i -x)fi/ii(x), 

with hi{x) analytic at x = uji. But then we deduce that (j) has all of the properties defining C , 
and since the closure of H{g) is symmetric on C' , it cannot have an eigenvector with nonreal 
eigenvalue. Thus = and the core property follows. 

After this digression we proceed to obtain our first main result, for which we again need to 
require g € 11,.. We denote the domain of the self-adjoint closure of H{g) by D(g), but keep the 
notation H(g) for the closure. We now use the key identity (jl.l7p as in the proof of Lemma [2.31 
to deduce 

{(t>,H{g)Iig)i;) = {H{g')I{g')(t),^), V,^,V^ G (2.11) 
Choosing A G C with 9A > 0, this entails in particular 

{(t),{H{g)+\)l{g)ij) = {{H{g') + \)l{g')^,i^), y^eD{g), y^P^D{g'). (2.12) 
Now for (j) G D{g) and tj: G D{g') there exist e, f £ TC such that 

<P = (Hig) + X)-'e, = {H[g') + A)" V, 

so ()2.12p implies 

{e,I{g){H{g') + \r^f) = {l{g'){H{g) + Xr'eJ), Ve,/ G H. 
Hence we obtain the relation 

Iig){Hig')+\)-' = iH{g)+X)-'l{g). 
From this we readily deduce 

[{H{g)+X)-\T{g)]=0, 

where 

oo 

T{g) = l{g)l{gr = ^ vl{g)en{g) ® e„(<7). (2.13) 

71=0 

The upshot is that the self-adjoint trace class operator T{g) commutes with the H{g)- 
resolvent. Since T[g) has trivial kernel by (j2.4p . all of its eigenspaces are finite-dimensional, 
and the //(^f) -resolvent leaves these finite-dimensional spaces invariant. Thus we can choose the 
eigenvectors of T{g) to be eigenvectors of the ff((7)-resolvent, hence of H{g) as well. We now 
summarize and extend these findings. 

Theorem 2.4. Let g £ H^- Then the positive trace class operator T{g) (|2.13p has an ONB 

{cn(5)}neN; of eigenvectors that are also eigenvectors of H[g). The functions en{g;x) are real- 
valued and belong to D{gQ,gi). They are of the form ()2.6p . where dn{g; x) extends to a function 
that is analytic in < a/2, even and n/r -periodic. Moreover, we have 

d„(<7;0)/0, dn{g;n/2r)^0, (2.14) 



and the spectrum of H{g) is nondegenerate. 
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Proof. We have just proved the first assertion. The second and third one have been shown 
earher too, cf. the paragraphs containing (|2.6|) and (j2.8|) . To prove the last one, we first note that 
it follows from en{g) S D{go,gi) that e„((7; x) is a classical solution to the ODE H(^g^en — -^n^n 
for a certain G M. From the Frobenius method it then follows that dn{g',x) satisfies (j2.14p . 
Moreover, any solution to the same ODE that is linearly independent of en{g;x) does not have 
the features encoded in (12. 6p . so the eigenvalue An is simple. ■ 

Quite likely, for G n\nj. the i/(5)-eigenvectors still belong to D{gQ,gi) and have the same 
features as for g G 11^. (If so, it follows as before that a[H[g)) is simple.) However, we cannot 
prove this via the above arguments, since they hinge on exploiting the HS operator T{g) with g 
restricted to 11^. 

It is a corollary of Theorem 12.41 that there exist distinct real numbers \n{g) such that 

H{g)en{g)=\n{g)en{g), riGN, g ^Tir- (2.15) 

At this stage, however, we have no other information on the -ff(5)-eigenvalues. In the following 
theorem, we obtain further results on the -ff((7)-spectrum that do not involve 1{g), and which are 
valid for all 5 G IT. As a preparation, note that H{g) is obviously bounded below for go,gi > 1. 
(Indeed, this entails that V{g; x) is bounded below, cf. (|1.2p .) For go and/or gi in (0, 1), however, 
V{g;x) diverges to —00 as x | and/or x | 7r/2r. Even so, our next theorem implies that H{g) 
is still bounded below. Its proof is based on a comparison argument involving the trigonometric 
BCi Calogero-Moser Hamiltonian (jl.2ip . 

Theorem 2.5. For all 5 G 11, the spectrum of the operator H[g) is a discrete set that is bounded 
below. Moreover, each eigenvalue has finite multiplicity. 

Proof. The trigonometric and elliptic differential operators (jl.2ip and (jl.ip differ by a real- 
valued potential 

yd{g]x) = Vt{go,gi;x) - V{g;x), 

that is bounded on [0,7r/2r]. Hence we can associate to ()1.2ip the operator on D{gQ,gi) given 
by 

Ht{go,gi) = H{g) + Vd{g;-). (2.16) 

This operator is essentially self-adjoint on D{gQ,gi), since H{g) is e.s.a. on D{gQ,gi) by virtue 
of Lemma 12.31 and the multiplication operator Vd{g; ■) is bounded and self-adjoint. We now 
introduce the trigonometric weight function 

wt{go, gi',x) = (sinrx)^^"(cosrx)^^\ 

and observe that we may rewrite D{gQ,gi) as 

D{go,gi) = wt{gQ,gi-xf/'^V. 

(Indeed, V (12. 8p is closed under products and we have 

s{x)\^ (s{./2r-x)Y ^^^^^^ 



sm rx I \ cos rx 



as is easily checked.) 

Next, we calculate the similarity-transformed differential operator 

At{gG,gi]x) = wt {go ,gi\x)~^/'^Ht {go ,gi;x)wt{go,gi; x) 
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= - 2r(5ocot rx - gitanrx)— + r'^{go + gif . 

This yields an essentially self-adjoint operator At{gQ,gi) on the dense subspace V (given by (|2.8p ) 
of the Hilbert space 

Ht = L^il^^TT /2r],wt{gQ, gi]x)dx). 

The crux is now that At {go , gi ) has an orthogonal base of eigenvectors belonging to P with 
nonnegative eigenvalues, explicitly given by 

dm,t{9o,9i;x) =2Fi{-m,go+gi+m;go + l/2;sm^rx), m G N, (2.17) 
At{go,9i)dm,t{9o,gi) = r'^{9o + 9i + 2mfdm,t{9o,9i), m G N. (2.18) 

Indeed, the r.h.s. of (j2.17p yields polynomials of degree m in sin^ rx, and these polynomials 
amount to the Jacobi polynomials 

Pt^Hy), a = -1/2, P = gi-l/2, y = cos2rx, m G N, 
which form an orthogonal base in 

np^L^[-l,l],{l-yni+yfdy), 
see for instance pO] . 

From the above it follows that the spectrum of Ht{go,gi), go,gi > —1/2, is bounded below, 
so by (j2.16p the same is true for H{g), g € H. Furthermore, it is clear from the eigenvalues 
(|2.18p that the resolvent of Ht{go,gi) is compact. By the second resolvent identity this also 
holds true for the resolvent of H{g), so the theorem follows. ■ 

This theorem also reveals the special character of the (at most four) self-adjoint extensions 
of H{g) restricted to C (j2.10p that arise via 1(g)- They are the extensions whose trigonomet- 
ric limits yield the (at most four) self-adjoint extensions of Ht{go,gi) restricted to C that are 
associated with the Jacobi polynomials. Of course, this is particularly clear when V{g; x) and 
Vt{go, gi] x) vanish identically, cf. Section HI 

For the remainder of this section we assume g G 11,. . It follows from Theorems 12.41 and 12.51 
that H{g) has nondegenerate eigenvalues 

Eo{g) < Ei{g) < E2{g) < ■ ■ ■ . 

Also, since H{g) and Ht{go,gi) differ by a bounded self-adjoint operator (cf. (|2.16p ). it follows 
from (I2.18P that the eigenvalues have asymptotics 

Emig) ~ Ar'^m'^, m — > oo. 

By Theorem 12.41 the associated eigenvectors Um{g) are also eigenvectors of T{g), and they are 
just the scalar multiples of en{g) for some n G N. Hence we may and will fix the ordering 
choice of the r((7)-eigenvectors for a given degenerate eigenvalue (if any occur), attained for 
example for i^j = Vj+i = • • • = z^j+fc, by requiring that the associated i?((7)-eigenvalues satisfy 
\j < Aj+i < • • • < \j+k ■ With this ordering of the en{g) understood from now on, there is 
a uniquely determined permutation Tg of the nonnegative integers such that 

K[m){9) = EM, 'T^eN, (2.19) 



cf. (I2T5D . 
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Next, we observe that we stih have not completely fixed the choice of the basis vectors en{g)- 
Indeed, we have required real-valuedness of en{g',x), but this leaves a sign undetermined. We 
may and will fix this ambiguity by requiring 

dn{g;0)>0, VnGN, 

cf. (EH]). 

At this point the alert reader might (or rather should) object that the sign choice just detailed 
may be in conflict with the decomposition (j2.5p . Indeed, once the ONB-functions en{g',x) have 
been fixed, the ONB-functions en{g'\x) are uniquely determined, and the signs of the latter 
might differ from the above sign convention. 

It follows from the last theorem of this section that this contingency does not arise. Its main 
result is, however, that H[g) and H{g') have the same spectrum. 

Theorem 2.6. The integral operator T[g) is real- analytic in g on li^ in the Hilbert-Schmidt 
norm. With the above sign convention for en{g;x) understood, it has singular value decomposi- 
tion (I1.19P . Furthermore, we have an equality of spectra 

a{H{g))=a{H{g')), geUr. 

Proof. Analyticity oiZ{g) in g with respect to the HS norm amounts to analyticity of its ker- 
nel ^{g;x,y) in g in the strong topology of the Hilbert space TCk (I1.14p . From the explicit for- 
mula (|1.4p it is readily verified that the kernel is strongly analytic in in a complex neighborhood 
of YLr, so the first assertion follows. 

Next, we fix g S 11^ and note that the vectors 

bm.{g;x) = e^^(^rn){g;x), m G N, 

yield an ONB with i7(5)-eigenvalues Em{g), cf. ([2T9]) . Recalling (pTT]) and (p3]) . we now 
deduce 

H{g')I{g')bm{g) =I{g')H{g)bm{g) = Em{g)T{g')bm{g). 

In view of (12. 4p . this entails that the vectors I{g')bm{g) are eigenvectors of H[g') with eigenvalues 
Em{g), and that they are total in 7i. Hence, 

Emig) = E„,ig'), Vm G N, 

and so the last assertion follows. Moreover, it follows that we have 

and 

Ag%n{g) = l^m{g')brn{g'), flm{g') G M*. 

To prove the theorem, it remains to show that all of the numbers fim{g) are positive. Now 
they do not vanish and are continuous on 11,., since their squares yield eigenvalues of the self- 
adjoint trace class operator T{g), which (by the first assertion) is continuous in g on 11,. in the 
trace norm topology. As 11^ is a connected set, it suffices to check positivity for one in 11^. For 
g = (0,0, 1, 1) we have g' = g and Tig) is self-adjoint. Thus the numbers ^m{g) are eigenvalues 
of 1{g)- These are easily determined explicitly, and they are all positive, cf. (j4.ip . ■ 

For all of the special g E 11,. where we hajVe ca,lculajteci the numbers 

(g) explicitly (cf. Sec- 
tion [4]), they are distinct and decreasing in m: 

Ma) > Mg) > Mg) > • • • > o. 

This entails in particular that the permutation Tg is the identity permutation. It is an open 
problem whether these features hold true on all of 11,. . 
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3 A hidden S4 symmetry 

In this section our purpose is to study, for a given coupling vector g € Hr, the couphngs g for 
which the spectra of H{g) and H{g) are equal. To this end we note that we can establish equality 
of spectra by exploiting the core property of D{gQ,gi), cf. Lemma 12.31 First, since the action 
of the differential operator H{g) on the domain D{gQ,gi) is invariant under taking gt ^ 1 — gt, 
and the domain does not depend on 52 and 53, we get equality of spectra (and even equality of 
operators) under the transformations 

52^1-52, 53^1-53- (3.1) 
(A caveat is in order: these maps can take g out of 11^.) Less trivially, consider the permutation 

(90,91,92,93) {91,90,93,92)- (3.2) 
Defining the 'mirror' unitary 

(M/)(x) = /(7r/2r-x), xe[0,7r/2r], f e H, (3.3) 
we clearly get 

MD{go,gi) = D{gi,go). 
This implies 

MH{go,gi,g2,g3)M = -H'(5'i, ffo, 53, 92), 

so that (j3.2p yields operators with equal spectra; in this case the two operators are different for 
generic g. 

Next, we note that the matrix Jjy (|1.5|) can be viewed as the reflection with respect to the 
unit vector (1, —1, 1, — 1)/2: 

Jat = I4 - 2'u (g) -u, -u = (1, -1, 1, -l)/2. 

Thus the self-dual couplings are given by 

9=9 ^ 90 + 92 = 91 +93, 

and on this 3-dimensional submanifold of 11^ the HS operator T[g) is self-adjoint. The maps (j3.ip 
and (j3.2p just considered do not involve Jn- But it is clear that even when we start from a self- 
dual g, the maps (|3.ip give rise to new couplings g that are generically not self-dual. When also 
^ E n^, we can act with to obtain g' E 11,. for which H[g') has the same spectrum as H{g) 
(by virtue of Theorem 12. 6p . 

All of the maps studied thus far are involutions. The question now arises which group is 
generated by these involutions. We proceed to answer this question. 

For this purpose it is expedient to switch to different couplings. We do this in two steps. 
First, we set 

A = g-C/2, C = (1,1, 1,1), (3.4) 
so that 

H{g) - + E - i) + (3-5) 
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Also, dMD and ([3:2]) become 

A2 -A2, A3 -A3, (3.6) 
(Ao, Ai, A2, A3) (Ai, Ao, A3, A2), (3.7) 

and we have 

A = A' 4^ Ao + A2 = Ai + A3, 

since C is self-dual. 

The second step may seem unmotivated at first sight. It consists in defining couplings 

co = Ao + A3, C3 = Ao-A3, ci = Ai + A2, C2 = Ai - A2. (3.8) 
The point of this definition is that the maps (j3.6p become permutations 

Co ^ C3, ci ^ C2, (3.9) 
while (j3.7p turns into the permutation 

(co, Ci, C2, C3) (ci, Co, C3, C2). (3.10) 
Moreover, the JAr-action transforms into the permutation 

C2 ^ C3, (3.11) 
as is easily checked. 

Since the three transpositions (j3.9p and (j3.1ip already generate S'4, it follows that (j3.10p 
may be viewed as a suitable product of these maps. The symmetries (j3.9p . (|3.10p generate 
an 8-element 54-subgroup isomorphic to the dihedral group /2(4). This subgroup gives rise to 
a decomposition of ^4 into three cosets. Specifically, as coset-representants we may choose the 
three permutations 

(co,Ci,C2,C3) (co,Ci,C2,C3), (co , Ci , C3 , C2) , (co , C2 , C3 , Ci ) . (3.12) 

Now that we have clarified the character of the symmetry group by exploiting the new cou- 
plings c, we may and will transform back to g. In particular, for a given g the three maps (|3.12p 
are 

9 9,g',g, 

where 

1/ ,1, 
90 = 2 (fo + 51 + 52 + 53 - 1), 91 = ^[90 + 91 -92 - 93 + 1), 

92 = 2^-90 + 91 -92 + 93 + 1), 93 = - 9i - 92 + 93 + !)• 

We now summarize and extend the above analysis in the following theorem. 
Theorem 3.1. The group G generated by the maps 

(50,51,52,53) (50,51, 1 - 52,53), (3.13) 
(50,51,52,53) (50,51,52,1 -53), (3.14) 
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and 



(3.15) 



with Jjsf given by (jl.Sp . is isomorphic to 5*4 and contains the mirror map 

(90,91,92,93) {91,90,93,92)- 

Now assume g belongs to the set IIg, cf. ()1.22p . Then the G-orbit of g belongs to Hg and one 
has 



Proof. We have already shown the first assertion. Next, we recall that 11^ is defined by the 
five linear constraints 

2go > -1, 2^1 > -1, 9o + 9i + 92 + 93 >0, 9o + 9i - 92 + 93 > 
9o + 9i+92 - 93 > -1- 

Together with the extra constraint defining IIg, i.e., 

90 + 91 - 92 - 93 > -2, 

this yields six constraints that transform into 

Cf, + c^>-2, G {0,1,2,3}, n^u, (3.17) 

as is readily verified from (j3.8p and (j3.4p . Since the latter are manifestly permutation invariant, 
the G-orbit assertion follows. Finally, since (|3.16p holds for the generators (|3.13p ~ (|3.15p and Hq 
is invariant under G, it holds for all w £ G. ■ 

Rewriting (|3.5p in terms of c, we obtain 



H{c) = --f^ - p{x + UJt) + J {cl + cl) p{x + ut) + -coC3{p{x) - p{x + ^3)) 



Requiring (j3.17p . the spectral invariance under the /2(4) subgroup is readily understood in this 
guise as well. But the invariance under any other c-permutation looks bizarre. 

4 Special cases 

It is plain from p.ip . (jl.2p that the Heun differential operator reduces to —d'^/dx'^ for the 
sixteen g's obtained by choosing gt £ {0, 1}, t = 0, 1, 2, 3. Viewing jdx^ as an operator on H 
with domain D{gQ,gi), it is also clear which ONB {6n}nGN of eigenvectors in D{gQ,gi) arises. 
Specifically, with the normalization constant 



^iH{g)) 



a{H{w{g))) 



\/w G G. 



(3.16) 



3 




AA= (4r/7r)i/2, 



the following ONBs are involved: 



(i) 50 = 51 = ^ 60 = AA/ V2, b, 



AAcos 2nrx, n > 0, 
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(ii) go = gi = I =^ 6„ = AAsin(2n + 2)rx, 

(iii) go = 1, gi = ^ bn = AAsm(2n + l)rx, 

(iv) go = 0, gi = 1 =^ bn = M cos{2n + l)rx. 

We proceed to study the above (7-choices in relation to 2{g)- For case (i) we must choose 

52 = 53 = 1 to obtain g £ 11^. (For 52 = 53 = we get Sg = 0, and for (72 = 0, (73 = 1 and g2 = 1, 

53 = we get g'l = — 1/2 and gQ = —1/2, resp.) Then we obtain from (jl.4p - (jl.6p and (|1.13p 

^-((0, 0, 1, 1); X, y) = R{x)R{uJi - x)R{y)R{uJi - y)/R{x + y)R{x - y). 

It now follows from Theorem 12.41 that the ONB (i) consists of eigenvectors for X((0, 0, 1, 1)). 
In fact, this can be readily verified by a contour integration, which also yields the eigenvalues 
explicitly. (To simplify the residues, one needs the duplication formula for the i?-function, 
cf. (jl.lOp .) Specifically, we obtain 

i/„ = 7r/pcosh(nra), g = (0,0, 1,1), c = (0, 0, -1, -1), (4.1) 

with p the infinite product in (|l.lip . 

Turning to case (ii), we get two self-dual subcases. For the first one, 

^((1, 1,0, 0); X, y) = /s(x)s(a;i - x)s{y)s{ui - y)/R{x + y)R{x - y), 

the eigenvalues of the corresponding integral operator I{g) on the ONB (ii) can be calculated 
in the same way as before. The result is 

i/„ = 7re™/pcosh((n + l)ra), 5 = (1,1,0,0), c= (0,0,1,1). 

Using (jl.lOp and (jl.lip . the second subcase can be rewritten as 

^r((l, 1, 1, 1); X, y) = ph{2x)s{2y)/R\x + y)R\x - y). 

With due effort (involving in particular a suitable use of (jl.lip ). the eigenvalues for this (7-choice 
can again be determined explicitly, yielding 

27r(n+l)re™ , , , , 

g = (l,l,l,l), c=(l,l,0,0). 



sinh((n + l)ra) ' 

It should be noted that the T((7)-spectra of these two subcases do not coincide, even though they 
are related by a c-permutation. 

The remaining two subcases are not self-dual, but they are related by mirror symmetry, 
cf. (j3.3p . Hence we only study 

5 = (1,1,0,1), c=(l,0,l,0) 4^ 5' = ^(3,1,1,1), c' = (1,0,0,1). (4.2) 

For this choice it is clear that H(g') has a nontrivial potential. Its eigenvector ONB is related 
to the ONB (ii) vial{g'). Specifically, we must have 

6„(5';x)=AA„ r^'^ dy [7^^^ ^^^'(3 sin(2n + 2)ry, g = {1,1,0,1), 

Jo [R[x + y)R{x-y)]-^/^ 

where the normalization constant Mn ensures that bn{g']x) has norm 1. Also, the H[g')- 
eigenvalues are given by 

E„((3,l,l,l)/2) = (2n + 2)2r2. (4.3) 
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The cases (iii) and (iv) are related by mirror symmetry, so we only consider case (iii). First 
we choose 

5 = (1,0,1,1), c=(l,0,-l,0) ^ 5' = ^(1,1,1,3), c' = (1,0,0,-1). 
Thus H(g') has a nontrivial potential. Just as for ()4.2p . its eigenvector ONB is given by 

6„(5';x)=AA„ r^'^ dy ^^^'(3 sin(2n + l)ry, 5 = (1,0,1,1). 

Jo [R{x + y)R{x-y)\y'^ 

The analog of (j4.3p is 

E„((l,l,l,3)/2) = {2n + l)W 
We continue with 

<7 = (1,0,0,0), c=(0,-l,0,l) ^ 5' = ^(1,1,-1,1), c' = (0,-1,1,0). 
Here we obtain the H{g')-ONB 

bn[g;x)=Mn dy— — ■ — — — ^ sm(2n + l)ry, 5 = (1,0,0,0), 

Jo [R{x + y)R{x-y}]^''^ 

with eigenvalues 

E„((l,l,-l,l)/2) = {2n + l)\^ 
Next we consider the self-dual choice 52 = 0, 53 = 1. Then we get 

^-((1, 0, 0, 1)52;, y) = ph{x)R{uJi - x)s{y)R{uJi - y)/R{x + y)R{x - y). 
As before, the eigenvalues of T{g) on the ONB (iii) can be explicitly determined, yielding 

iy„ = 7re™/Vpcosh((ri + l/2)ra), 5 = (1, 0, 0, 1), c = (1, -1, 0, 0). 
Finally, we study the choice 

<7 = (1,0,1,0), c = (0,0,-l,l) ^ ^' = (0,1,0,1), c' = (0,0, 1,-1), 
for which 

^-((1, 0, 1, 0); X, y) = ph{x)R{x)s{oJi - y)R{ui - y)/R{x + y)R{x - y). 
Once more, a contour integration yields the expected result 

/>7r/2r 

sin(2n + l)rx = f„ / (iy^((l, 0, 1, 0); x, y) cos(2n + l)ry, 
Jo 

with the singular values Vn given by 

zy„ = 7re™/Vpsinh((n + l/2)ra), 9 = (1, 0, 1, 0), c = (0, 0, -1, 1). 

We conclude this section with some remarks on special cases of a different character. We 
recall first that I{g) reduces to a rank-one operator for 5 G 11 with Sg = 0, cf. p.4|) and (|1.6|) . By 
continuity in g, it therefore follows that the functions w{g; x)^/^ and 'w{g';x)^^'^ are eigenfunctions 
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of H{g) and H{g'), resp., with the same eigenvalue. We point out that this is not at all obvious. 
With due effort, however, this can be directly verified. 

Specifically, it follows from I (3.79), (3.113) and (3.114) that the eigenfunction property is 
equivalent to constancy of the function 

We encountered this functional identity before on p. 253 of [18], where we supplied a (short) 
proof. Using next I (3.105)-(3.106), the constant (which amounts to the ground state eigenvalue) 
can be calculated explicitly: 

3 3 

Eo{g) =^gjej{gj + 2gQ), ej = p{ujj), j = 1,2,3, ^gt = 0. 
j=i t=Q 

From this formula one readily verifies equality to EQ{g'). 

Finally, we comment on the case where V{g; x) is real-analytic on M. For this we need to 
choose 50,51 S {0,1}, cf. (11. 2p . Then one need only combine the W-ONB's bn{{0,0, g2, gs); x) 
and 6n((l, 1, 52, 53); a;) (consisting of real-analytic vr/r-periodic functions that are even and odd, 
resp.) to obtain an ONB of vr/r-periodic functions for 

Tib = L'^i[-TT/2r,Tr/2r],dx). 

Likewise, an ONB of vr/r- ant iper iodic functions for Tib results from combining the "H-ONB 
6„((0, 1, 52, 53); a^) of even functions and the "H-ONB of odd functions 6„((1, 0, 52, 53); 2;). In 
terms of Floquet theory, therefore, one can only arrive at the Hfo-ONB's associated with the 
multipliers 1 and —1. Whenever the four 5- values are all in 11^, one can invoke the spectral 
invariance under taking g g' . Note that V{g'; x) may be singular at x = and/or x = 7r/2r, 
as we have already seen in the above special cases where V{g; x) vanishes identically. 

A Proofs of Lemmas 12.11 and 12.21 

Proof of Lemma 12.11 We recall that the function 

/•■K/2r 

h{x)= dyS{g;x,y)^{y), (j) ^ L\%Ti /2r\), xe[0,^/2r], 







extends to an even analytic function in llmzj < a/2, cf. the paragraph containing (12. ip . Thus 
its vanishing is equivalent to vanishing of the coefficients in its Taylor expansion 

00 

h{x) =Y,Ckx''\ (A.l) 

A:=0 

at x = 0. We proceed to study these coefficients. 

Obviously, we need information on the coefficients in the expansion of S{g;x,y) at x = 0. 
To this end we first consider the expansion 



^ 00 
F{x) = - - (In R{x + y) + In R{x - y)) = ^ 



(2j)! 
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Using 



"dz^ l'^-^(^) = P(^) + 2r/r/7r, p{z) = p{z + ia/2), 



(cf. (fTTTI) . (frT2l) ). we obtain 

F(0) = -lnii(y), F(2)(0) = p{y) + 2r]r / n , 
F(2'=)(0) = p(2A:-2)(^)^ A;>1. 

Now from standard elliptic function lore [1] we deduce 

p(2fc-2)(^) = (2A; - iy.p{zf + l.d., > 1, 

where l.d. stands for a polynomial in p{z) of degree < k. Thus we have 

oo 

F{x) = - In R{y) + ^ a^x^^ 

with 

ak = ^p{y)'' + Id. (A.2) 
We are now prepared to consider the expansion of 5. It is given by 

(oo \ oo 

j=l J k=0 

so that the first few bk read 

6o = 1, ^1 = 2sgai, 62 = 2sga2 + 2Sgai, 63 = 25^03 + 45^0201 + iSgai/S. 

More generally, it is easily seen that bk is the solution to the recurrence 
k 

kbk = 2sg'Yjajbk^j, k > 0, bo = I. 

In view of (IA.2|) . this implies that bk is given by 

bk = Q'-'Ksg)p{y)'' + Id., (A.4) 
where Q^^\u) is a polynomial of the form 

k 



1=1 

We are now in the position to use our assumption Sg > 0. It entails that the coeffi- 
cient Q^'^^Sg) in ()A.4p is positive: 

bk = qkp{y)^ + l-d., gfcG(0,oo), V/cEN. (A.5) 

Returning to (jA.ip . we first note that ()A.3P entails 

/>7r/2r 

Cfc = / dybkiy)ipiy), 
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where 

V^(y) ^ R{y)-^'^(p{y). (A.6) 
Since the coefficients vanish by assumption, it follows recursively from ()A.5p that we have 

r>7r/2r 

dyp(y)V(y) = 0, VfcGN. 





Now p{y) is monotone increasing on [0,7r/2r]. Thus we can invoke the Stone- Weierstrass 
theorem [15] to conclude that the span of the functions 1, p(y), p(y)^, p(y)^, ... is dense (in the 
supremum norm) in the space of real- valued continuous functions on [0,7r/2r]. Since ip{y) is an 
L^-function, it follows by dominated convergence that we have 

■k/2t 

dyC{y)iP{y) = 0, VC G Cm([0, 7r/2r]). 





Hence ^ = 0, so by (IXB]) (/> = 0. ■ 

Proof of Lemma 12.21 Letting e{x) G D{gQ,gi), we first show that the action of the differential 
operator H{g) (II. ip on e{x) yields a function in 7i. To this end we use the relations 

eix) = sixy°s{x-u;iy'd{x), d£V, 

and 

s"{z) 



s{zy s{z) 

to calculate 



t=0,l 



, , ( ^, , dl{x\ s'ix^uji^ d"ix) 



where 



Q{x) = 2goffi ; ! ; (a.?) 

S[X) S[X + LOl) 

Subtracting the potential term 

3 

e{x)'Ygt{gt - l)p{x + ut), 
t=o 

we obtain a sum of terms that are manifestly in 7i, save for the terms 

e(x) (q{x) + Y: I + 2K-)^°Kx - u^y^ Y 9td'{x)'-^j^±^. (A.8) 

Now s(2;) is odd and has a first order zero for 2; = 0, so s"{z) vanishes at z = at least to 
first order. Thus the two terms in the first sum are smooth for x G M. Next, s'{z) vanishes at 
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z = iOi, since s{z) is odd and 2u;i-antiperiodic. Therefore Q{x) ()A.7p is also smooth. Finally, 
since d'{0) = d'{uji) = 0, the two terms in the second sum are continuous on [0,Ci;i]. Hence the 
function given by (jA.SP is in 7i, so that H{g) is well defined on D{gQ,gi). 

It remains to show that H{g) is symmetric on D{gQ,gi). To this end, let di,d2 £ T) and 
consider 

- {H{g)s{-f<^s{- + s(.)5°s(- + ^iTd2). (A.9) 

From the above calculation we see that this equals 







dxs{xf30s(^rj. + ^^fgi (x) - di {x)d2 (x) 

+ 2[(i';L(3^)'^2(x) - di{x)d'2{x)] ^ 5r'^ '^^ ~^ '^'^ 



Clearly, this can be rewritten as 
r'dx (-d^-^[s{xf30s^ 

X + LOl f9^d'^{x)]+d2{x)-^[s{xf30s{ X + LOi)'^^'^ d[{x)] I . (A. 10) 
Jo \ dx dx J 

Now for all d £ D we have 

d'{x) = / d"{t)dt = 0{x), X I 0, d'{x) = - d"{t)dt = 0{x - wi), x t wi- 



Integrating by parts in (jA.lOp . it follows that the boundary terms vanish, whereas the new 
integrand vanishes identically. As a result, the difference (IA.9j) vanishes, so that H{g) is indeed 
symmetric on D (go, B 
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